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Abstract 

We utilise the method of Wiener reconstruction with spherical harmonics and Bessel functions to recover 
the density and velocity fields from the IRAS 1.2Jy redshift survey. The reconstruction relies on prior 
knowledge of the IRAS power-spectrum and the combination of density and bias parameters f3 = il - 6 /b. 
The results are robust to changes in these prior parameters and the number of expansion coefficients. Maps 
of these fields are presented in a variety of projections. Many known structures are observed, including clear 
confirmations of the clusters N1600 and A3627. The Perseus-Pisces supercluster appears to extend out to 
~ 9000 km s _1 , and the reconstruction shows 'backside infall' to the Centaurus/Great Attractor region. A 
qualitative comparison of the reconstructed IRAS gravity field with that from Tully-Fisher peculiar velocity 
measurements (Mark III) shows reasonable agreement. The Wiener reconstruction of the density field is also 
the optimal reconstruction (in the minimum variance sense) of any quantity which is linear in the density 
contrast. We show reconstructions of three such quantities. The misalignment angle between the IRAS and 
CMB Local Group dipoles is only 13° out to 5000 km s~\ but increases to 25° out to 20,000 km s _1 . The 
reconstructed IRAS bulk flow out to 5000 km s _1 is ~ 300 km s _1 , which agrees in amplitude with that 
derived from the Mark III peculiar velocities (~ 370 km s" 1 ). However, the two bulk flow vectors deviate 
by ~ 70°. Finally, moment of inertia analysis shows that the Wiener reconstructed Supergalactic Plane is 
aligned within ~ 30° of that defined by de Vaucouleurs. 

1. Introduction 



Observations of large scale structure have been advancing spectacularly. Both in terms of redshifts and 
galaxy peculiar velocities, the number of measurements has grown by more than a factor of ten over the 
past decade. Hand in hand with observational progress, a plethora of statistical methods has been developed 
to analyse the data, each attempting to extract as much astrophysical information as possible. Accurate 
determination of the galaxy distribution and velocity field would not only place constraints on cosmological 
parameters, but would also provide an insight into the mechanisms of structure formation which generated 
the complex pattern of sheets and filaments we observe today. 

Redshifts have been determined for many more galaxies than those for which we have direct distance 
measurements. Over the next few years, automated surveys such as the 2-degree-Field (2dF) and Sloan 
Digital Sky Survey (SDSS) will increase this lead even further, as over a million new galaxy redshifts become 
available. This wealth of data has inspired a great deal of work on techniques for reconstructing the velocity 
field on the basis of redshifts alone. Given the assumptions of linear mass-to-light biasing, and purely linear 
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structure evolution, the mapping between real and redshift-spacc is well denned (Kaiser 1987). However, 
in applying reconstruction methods there is freedom of choice regarding (i) the functional representation 
(e.g. Cartesian, Fourier, Spherical Harmonics, or Wavelets) and (ii) the filtering, or smoothing scheme (e.g. 
a Gaussian sphere, a sharp cut-off in Fourier space, or a Wiener filter). Smoothing is necessary to reduce 
sampling noise before the real-space density field can be calculated from that in redshift-space. 

The first reconstruction technique to be applied to redshift surveys was based on iterative solution of the 
equations of linear theory, pioneered by Yahil et al. (1991), hereafter YSDH (see also Yahil 1988; Strauss 
& Davis 1988). This involves solving for the gravity field in redshift space, which can then be used to 
derive an estimate of the peculiar velocities for a given value of f3 = £l 06 /b. These velocities allow the 
rcdshifts to be corrected, providing an updated set of distance estimates. This is repeated until the distance 
estimates converge. YSDH use variable smoothing, but their smoothing regime is not rigorously formulated. 
Variants of this technique have been successfully applied to IRAS selected galaxy catalogues (YSDH; Kaiser 
et al. 1991). 

The reconstruction procedure used in this paper follows Fisher et al. (1995a; hereafter FLHLZ) and is 
non-iterative. The density field in redshift-space is expanded in terms of spherical harmonics and Bessel 
functions. The real-space density, velocity and potential fields are reconstructed using linear theory and a 
Wiener filter which assumes a given power spectrum and noise level. It is important to note that, as opposed 
to ad-hoc smoothing schemes, the smoothing due to a Wiener filter is determined by the sparseness of data 
relative to the expected signal. 

Cosmography using spherical harmonics is not a new technique (Peebles 1973), but it has enjoyed a 
renaissance with the advent of near whole-sky galaxy surveys (e.g. Fabbri & Natale 1989; Regos & Szalay 
1989; Lynden-Bell 1991; Scharf et al. 1992; Scharf & Lahav 1993; Fisher, Scharf & Lahav 1994, hereafter FSL; 
FLHLZ; Lahav 1994; Nusser & Davis 1994; Heavens & Taylor 1995). Similarly, the Wiener filter (Wiener 
1949; hereafter WF) has long been used in engineering to recover the best estimate of the true signal 
from one corrupted by imperfect measurement. Recently, Wiener filtering has been applied to a number 
of cosmological reconstruction problems. Lahav et al. (1994; hereafter LFHSZ) reconstructed the angular 
distribution of IRAS galaxies, while Bunn et al. (1994) applied the method to the temperature fluctuations 
in the COBE DMR maps. Kaiser & Stebbins (1991) used a similar method to reconstruct velocity fields. 
Reviews of the WF and linear estimation methods can be found in Press et al. (1992), Rybicki & Press 
(1992) and Zaroubi et al. (1995; hereafter ZHFL). 

FLHLZ applied the 3-dimcnsional Wiener reconstruction to the IRAS 1.2Jy survey, showing preliminary 
predictions for the density, velocity and potential fields. In this paper, we expand on the initial results in 
FLHLZ, presenting detailed maps of the reconstructed fields, as well as optimal determinations of the Local 
Group dipole, bulk flows and the extent of the Supergalactic Plane. We compare the IRAS reconstructed 
peculiar velocity field to that derived from Tully-Fisher distances (Mark III) and examine the sensitivity of 
the WF method to choice of parameters. 

In §2, we give a brief summary of the decomposition and filtering techniques used in this reconstruction 
and the data set to which they have been applied. §3 discusses the role of the prior model in the WF and 
the methodology behind selection of parameters. Results of the reconstruction are shown in §4, and these 
are compared to Mark III data in §5. As detailed below, the choice of prior model is central to the operation 
of the WF and §6 discusses the qualitative influence of these parameters on the reconstructed density and 
velocity fields. In §7 and §8, we show reconstructions of the Local Group dipole and bulk flow respectively; 
both determined directly from harmonic coefficients. Similarly, a reconstruction of the Supergalactic Plane 
is shown in §9, and its alignment compared to previous measurements. Finally §10 contains a discussion of 
the work in this paper. 



2. Wiener reconstruction of IRAS data 
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2.1 Reconstruction Method 



This paper concentrates on cosmography, presenting results from the application of the WF method to 
IRAS 1.2Jy survey data. Hence, we include here only a brief introduction to our reconstruction technique. 
The approach we use is discussed in detail by FLHLZ, while the formalism of Wiener filtering as it pertains 
to LSS reconstruction can be found in ZHFL. 

A density field, p, can be expanded in terms of spherical harmonics, Y| m , and Bessel functions, jf. 

(max +1 «max(0 

p(s) = ^2 E E Cln Plmn 3l{k n S)Y lm (s) . (1) 

Z=0 m— — l n—1 

The discrete k n 's are chosen according to the boundary conditions, so as to make the set orthogonal (sec 
Appendix A in FLHLZ). This process is analogous to Fourier decomposition, but instead using a set of 
spherical basis functions. The data from a rcdshift catalogue can be seen as a set of N discrete points, s,, 
each giving the direction and redshift of a galaxy. These are used to estimate the underlying density field in 
rcdshift space, p (s), expanded as in equation 1. Here, C;„ are normalisation constants, while the harmonic 
coefficients are given by 

N 1 

Pfmn = Y !J7^3i{Ksi)Y lm {s{) , (2) 

where <j>{si) is the spherical selection function of the survey, evaluated at the radius of the i th galaxy. From 
/5(s), we can derive an estimate of the overdensity field, 5(s) = (p(s)/p) — 1, in terms of the same basis 
functions, with coefficients of mn . This is not only more useful for dynamics, but also has, by construction, a 
mean of zero, which is useful later on in the reconstruction method. 

As redshifts are line-of-sight measurements, peculiar velocities will introduce only radial distortion in 
redshift space. Given that 6 s is expressed in terms of orthogonal radial and angular components, redshift 
distortion can be seen as coupling between the radial modes in real and redshift space, leaving the angular 
modes unaffected. This can be expressed by introducing a matrix which links radial modes in the density 
fluctuation fields, 

dfmn = (^)nn' ^hnn' > (3) 

n' 

where Z; depends on the chosen value of [3 and the selection function for the survey (see Appendix D in 
FLHLZ). Here we use superscripts R and S to denote real and redshift space respectively. Note that this 
is only applicable for all-sky surveys; in samples with incomplete sky coverage, there will also be more 
complicated coupling of angular modes. 

In a perfect galaxy catalogue, with arbitrarily high sampling, the real-space harmonics could be estimated 
by simple inversion of the coupling matrix in equation 3. However, with real data, straightforward inversion 
is often unstable, and a regularisation scheme is required. Without such a scheme, shot noise from the 
discreteness of the galaxy distribution might be amplified by inversion, leading to an estimate that was far 
from optimal. Application of a Wiener filter to the inversion process allows one to reconstruct the minimum 
variance solution, assuming a prior power spectrum and a shot noise level. Appendix A provides a summary 
of the WF method. For further details, please refer to Rybicki & Press (1992), LFHSZ and ZHFL. 

In our problem, the redshift coupling matrix represents the response function of the system, acting on 
the real-space distribution plus shot noise to give the observed, 5 s . Hence, the WF estimation of the real 
space overdensity field is 

(Owf = E ( S < [ S ' + ^"'L, ( Z T\> n » Sfmn" ■ (<0 

n'n" 

The assumed true signal matrix, S, depends on the power spectrum, while the noise matrix, N, is a function 
of the mean density and selection function (see LFHSZ). Incidentally, as Rybicki & Press (1992) point out, 
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the WF is in general a biased estimator of the mean field unless the field has zero mean. As mentioned 
above, the density fluctuation field has zero mean. 



Just as the density field can be decomposed into harmonics, so too can the radial and transverse com- 
ponents of the velocity field, and the potential field. In the case of radial velocities, linear theory allows 
the velocity due to inhomogeneities within r < R to be expressed directly in terms of the harmonics of the 
overdensity field. Again, this can be expressed by introducing a coupling matrix, in a similar manner to that 
in equation 3, such that 



Transverse velocity and potential harmonics are also related to those of the density field (Appendix C in 
FLHLZ). 

To summarise, the redshift-space density field is estimated in terms of orthogonal functions which are 
eigenfunctions of the Laplacian operator. For a given choice of /3, the radial coupling matrix linking real- and 
redshift-space harmonics can be determined. Rcgularisation by a Wiener filter, with a given power spectrum 
as its prior, allows inversion of the coupling matrix and determination of the real-space density field. Linear 
theory then allows the velocity and potential fields to be determined from the estimated density field. We 
emphasise that the method is non-iterative, and uses a filter which is optimal in terms of assumed signal and 
noise characteristics. It provides a non-parametric and minimum variance estimates of the density, velocity, 
and potential fields which are related by simple linear transformations. 

2.2 The IRAS 1.2 Jy Survey 

The IRAS surveys are uniform and complete down to Galactic latitudes as low as ±5° from the Galactic 
plane. This makes them ideal for estimating whole-sky density and velocity fields. Here, we use the 1.2 Jy 
IRAS survey (Fisher et al. 19956), consisting of 5313 galaxies, covering 87.6% of sky with the incomplete 
regions being dominated by the 8.7% of the sky with |6| < 5°. 

In principle, the method can be extended to explicitly account for the incomplete sky coverage. We 
adopt the simpler approach of smoothly interpolating the redshift distribution over the missing areas using 
the method described in YSDH. Given this, we can assume a purely radial selection function. In LFHSZ, we 
examined the validity of this interpolation and found its effect on the computed harmonics to be negligible 
for / £ 15 for the geometry of the 1.2 Jy survey. We therefore use the interpolated catalogue (and a simplified 
formalism) for our reconstructions. 

3. Selection of a prior model 

Throughout this paper, the density field is expanded within a spherical volume of radius R = 20, 000 km s _1 .| 
Boundary conditions are imposed at the edge of this volume, forcing the radial modes to be discrete rather 
than continuous. We have demanded that the logarithmic derivative of the gravitational potential be con- 
tinuous at the boundary; a discussion of various boundary conditions is contained in Appendix A of FLHLZ. 
We then truncate the summation over angular and radial modes at some Z max = 15 and k n R < 100 respec- 
tively, effectively determining the resolution of the density field (please see Appendix B of FLHLZ for further 
details). The qualitative effects of Z max on resolution are discussed in §6. 

Beyond the definition of the expansion, the reconstruction method itself has a number of free parameters. 
The WF assumes a prior power spectrum, while the redshift coupling matrix depends on choice of (3. 
Fortunately, the shape and amplitude of the power spectrum has been relatively well determined for IRAS 
galaxies. The power spectrum is described on scales ^ 200 /i _1 Mpcby a CDM power spectrum with shape 
parameter T (see Efstathiou, Bond, & White 1992) in the range r ~ 0.2 - 0.3 (Fisher et al. 1993; Feldman, 
Kaiser, & Peacock 1994). The normalisation of the power spectrum is conventionally specified by the variance 
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of the galaxy counts in spheres of 8 /i _1 Mpc, erg. Fisher et al. (1994a) used the projection of the redshift 
space correlation function to deduce the value, er 8 = 0.69 ± 0.04 in real space. We have adopted the WF 
prior given by T — 0.2 and as = 0.7. 

For simplicity, we assume linear, scale-independent biasing, where b measures the ratio between fluctua- 
tions in the galaxy and mass distribution, 

(<Wp)gai = K$P/P) 

mass (6) 

In linear theory, the density and velocity fields are related by f3 = fl 6 /b, where fi is the density parameter. 
The correct value of (3 for IRAS is uncertain. Strauss & Willick (1995) review work on determining /?iras 
by a variety of methods, giving values in the range 0.35 < /3iras < 1-28. 

As discussed in FLHLZ, our reconstruction technique is essentially pcrturbative, involving a Taylor 
expansion to first order in Av = (v(r) — v(0)), where v(0) is the motion of the observer at the same frame 
as v(r). Given that the choice of calculation frame is arbitrary, it seems sensible to select a frame in which 
Av is minimal. It is known that nearby galaxies (nearer than ~ 2000 km s -1 ) have small peculiar velocities 
relative to the Local Group (LG), while very distant galaxies are, on average, most likely to be at rest with 
respect to the Cosmic Microwave Background (CMB). Rather than taking an ad-hoc value for /3, it makes 
sense to try the reconstruction in both frames, and take the value of (3 for which the two calculations are in 
best agreement. 

We compare the two reconstructed line-of-sight peculiar velocities at the position of all catalogued IRAS 
galaxies by calculating the RMS scatter as function of (3. Figure 1 shows that the minimum scatter occurs 
at (3 ~ 0.7; this result is robust to changes in the limiting radius of the reconstruction. There are many 
other possible methods for selecting the prior value of (3. For example, we can choose the value which best 
reconstructs the motion of the LG with respect to the CMB; this gives (3 « 0.6 to recover the 600 km s _1 . 
Another alternative is to find the (3 for which the Local Group motion is identical in both calculation frame, 
yielding (3 ~ 0.8. FLHLZ used (3 = 1.0, based on maximum likelihood analysis in FSL. Hereafter we adopt 
f3 = 0.7 based on figure 1, which is well within the range quoted by various studies (Strauss & Willick 1995). 
In §6 we investigate the effects on the reconstruction of changing (3, T and a%. 

4. Density and Velocity Maps 

The maps shown in this section were derived by the technique detailed above. In summary, the recon- 
struction was carried out within a sphere of radius R = 20,000 km s , with angular modes limited to 
l max = 15 and radial modes to k n R < 100, with a total of 6905 coefficients. A canonical set of values were 
adopted for cosmological parameters, such that (3 — 0.7, T = 0.2 and a s = 0.7. The reconstruction has a 
variable smoothing scale similar to a Gaussian window of width, a s , proportional to the mean inter-particle 
separation (e.g. a s — 436 km s _1 at r — 4000 km s _1 , a s = 626 km s~ x at 6000 km s _1 , and a s = 1130 
km s _1 at 10,000 km s _1 ). The density contrast values shown in the diagrams and quoted below have been 
obtained under this smoothing regime. Figures 2 and 3 show Aitoff projections of the real-space density and 
radial velocity fields respectively, plotted in Galactic coordinates, evaluated across shells at various distances. 
Many previous papers have examined parts of the volume shown in these maps, identifying clusters and voids 
(e.g. Tully 1987, Pellegrini et al. 1990, Saunders et al. 1991). Where obvious analogues exist, structures are 
labelled in accordance with previous papers. 

Figure 2a shows a density shell at r = 2000 km s _1 . All major features have already been identified, 
and there is the expected concentration of clusters in the northern Galactic hemisphere. In particular, note 
Fornax-Doradus-Eridanus {I ~ 180°, b ~ -60°), N5846 (I ~ 350°, b ~ 45°), Virgo (I ~ 290°, b ~ 70°) and 
Ursa Major (I ~ 135°, 6-30°). The foregrounds of Hydra (I ~ 270°, b ~ 10°), and Centaurus (I ~ 310°, 
b ~ 20°) are visible, and these structures can be seen to extend onto the next shell. The observed shoulder 
on Hydra at / ~ 240°, b ~ 0° is due to Puppis (e.g. Lahav et al. 1993), which appears as a distinct cluster at 
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r = 2500 km s _1 . The ovcrdcnsity marked Ca (Z ~ 10°, b ~ —60°) has previously been identified as part of 
the Pavonis-Indus-Telescopium (P-I-T) supercluster (Santiago et al. 1995). However, in this reconstruction, 
Ca appears to be quite distinct, peaking at r ~ 2500 km s _1 with <5 max = 2.8 on the smoothing scale 
described above. The Local Void (LV), as discussed by Tully (1987) is clearly seen in the generally under- 
dense region roughly spanning 330° < I < 120°, —50° < b < 30°. Two additional voids are also marked, Va 
(150° < I < 220°, -30° < b < 30°) and V/3 (240° < I < 300°, -90° < b < -20°). The LV and V/3 regions 
extend out to r ~ 4500 km s _1 , while Va ends at r ~ 3500 km s _1 . 

Figure 2b is density at r = 4000 km s _1 . Continuing from the previous shell, the Hydra (Z ~ 270°, 
6 ~ 25°) and Centaurus/Great Attractor (GA; Z ~ 315°, b ~ 15°) superclusters are shown clearly, both 
extending across -10° < b < 45°. The foregrounds of P-I-T (Z ~ 340°, b ~ -25°) and Perseus-Pisces 
(P-P; Z ~ 150°, 6 ~ -15°) can be seen, as can Leo (Z ~ 240°, b ~ 75°), Cancer (Z ~ 190°, 6 ~ 25°) 
and Cetus (Z ~ 180°, 6 ~ -60°). There is also confirmation of N1600 (Z ~ 200, b ~ -28), which is 
the strongest overdensity on this shell (<5 max = 4.2). The position of N1600 closely matches that previously 
determined in both the IRAS (Saunders et al. 1991) and ORS (Santiago et al. 1995) surveys. Camelopardalis 
(I ~ 145°, b ~ 30°) is visible in the northern foreground of the Perseus-Pisces supercluster (centre at 
I ~ 150°, 6 ~ -15°). The two are linked by a cluster marked C/3 (Z ~ 135°, 6 ~ 10°), making a wall of 
continuous overdensity extending over —85° < b < 45°. Extensions of LV and Va are also seen, although 
the underdensity is considerably less smooth than at r = 2000 km s . Three dense regions on the shell 
do not correspond well to any previous labelled clusters. C7 (I ~ 110, b ~ —70) is the foreground of a 
structure which extends out to r ~ 8000 km s" 1 , including the previously-identified A194 galaxy cluster at 
<~ 6000 km s _1 . C(5 (Z ~ 105, b ~ —40) is the foreground of a cluster which extends to r ~ 6500 km s _1 , 
where it merges into Pegasus. Finally, Ce (Z ~ 65, b ~ 70) is the foreground of a weakly overdense wall that 
merges into the Great Wall structure at r ~ 8000 km s _1 . 

Figure 2c is density at r = 6000 km s _1 . Continuing from the previous shell, note clear extensions of 
P-P (Z ~ 150°, b ~ -15°), P-I-T (Z ~ 0°, 6 ~ -30°), Hydra (Z ~ 280°, 6 ~ 10°), Leo (Z ~ 270°, b ~ 60°), C 7 
(Z ~ 120°, 6 ~ -70°) and Ce (Z ~ 0°, b ~ 70°). Faint backgrounds can also be seen for Cancer (Z ~ 190°, 
b ~ 30°) and C5 (Z ~ 110°, 6 ~ -45°). Orion (Z ~ 190°, b ~ 0°), as shown on this shell, is the background of 
a strong cluster centred at r ~ 5000 km s , where it links Cancer and N1600 to form a structure extending 
over —45° < b < 45°. A569 (Z ~ 165°, b ~ 20°) is centred on this shell, and corresponds very well to 
the position, extent and redshift (Abell et al. 1989). The foreground of the Pegasus (Z ~ 90°, b ~ —20°) 
supercluster is seen weakly; in this reconstruction, Pegasus extends to r > 12,000 km s _1 . There is also 
confirmation of the cluster reported by Kraan-Korteweg et al. (1996), marked A3627 (Z ~ 120°, b ~ —8°), 
which can be seen distinctly now that Centaurus no longer swamps the region. The remaining cluster, 
labelled C( (Z ~ 40°, b ~ —20°) does not correspond well to any known labels; it is the foreground of a 
large structure extending to r > 10,000 km s _1 , and appears to be quite separate from the nearby P-I-T 
supercluster. 

Figure 2d is density at r = 8000 km s _1 . Continuing from the previous shell, P-P (Z ~ 150°, b ~ 0°), Leo 
(Z ~ 260°, b ~ 45°), A3627 (Z ~ 320°, b ~ -8°) and C( (Z ~ 30°, b ~ -35°) extend out tor- 8000 km s" 1 . 
The cluster marked Ce has merged into the body of the Great Wall (GW; (Z ~ 0°, 0° < b< 90°), which itself 
merges with C£ at r ~ 10, 000 km s _1 to form a wall extending over —60° < b < 90°. Another wall is formed 
by Leo and Coma (60° < Z < 240°, b ~ 75°); this also seems to extend down beyond the Galactic equator 
at greater distances. A779 (Z ~ 170°, b ~ 45°), A539 (Z ~ 190°, 6 ~ -20°) and A400 (Z ~ 170°, b ~ -40°) 
all correspond to clusters in the Abell catalogue. Cygnus (Z ~ 80°, b ~ 15°) and Cry (Z ~ 250°, 6 ~ ^5°) 
are both foregrounds of structures extending to r ~ 10,000 km s _1 . Notably, in this reconstruction, P-P 
extends out to roughly 9000 km s _1 ; far further than usually thought. Out at r > 6000 km s _1 , the peak 
overdensity for P-P lies very close to the Galactic equator, and well within the zone of avoidance (ZOA). It 
is important to remember that Cry and P-P both lie largely within the \b\ < 5° region in which data has been 
interpolated to fill the zone of avoidance. However, both are strong overdensities on this shell (c>p_p ~ 4), 
and the shoulders are clearly visible in the IRAS data. 

Figure 3a shows line-of-sight velocity at r = 2000 km s _1 . The major feature in this plot is the velocity 
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dipole due to the Centaurus/GA region (towards I ~ 300°, b ~ 15°). Also visible are the effects of the 
P-P/Cancer/N1600 mass concentration, causing outflow towards I ~ 165°, & ~ —10°. In figure 3b, at 
r = 4000 km s _1 , some backside infall is visible in the Centaurus/GA region. The strongest features are 
outflow towards P-P (I ~ 145°, b ~ -20°), N1600 (Z ~ 195°, 6 ~ -20°) and Cancer (/ ~ 190°, b ~ 20°), 
while weaker outflow can also be seem towards the centre of P-I-T (I ~ 340°, & ~ —20°). Figure 3c, at 
r = 6000 km s -1 is dominated by motion from the south Galactic pole towards the North Galactic pole. 
This is towards the Great Wall in the northern Galactic hemisphere, while outflow towards the A569/P-P 
region (I ~ 165°, —50° < b < 30°) and the overdensity around G-q (I ~ 240°, b ~ —15°) can also be seen. 
Backside infall is now seen behind P-I-T and southern Centaurus. Finally, in figure 3d, at r — 8000 km s _1 , 
we can see evidence of the mass distribution at larger radii. Beyond 10,000km s , this reconstruction places 
large superclusters at I ~ 60°, b ~ 60° (forming out of the Great Wall), and I ~ 220°, b ~ —30° (forming 
behind C?7 and A400). These cause the observed outflow, while strong infall can be seen behind Pegasus 
(I ~ 90°, b ~ -25°) and Hydra (I ~ 285°, & ~ 5°). 

Figure 4 shows the density field on planes at Supergalactic x, y, z — 0, ±4000 km s _1 , while figure 5 
shows the corresponding velocity field. The Supergalactic Plane can be clearly seen on the SGZ = km s _1 
slice, with Centaurus and P-P visible as strong overdensities. The velocity field on the same plane shows 
the tug-of-war between these two superclusters, and the resultant effect on the local group. There has 
been some controversy over the presence of backside infall in the GA region (e.g. Dressier & Faber 1990, 
Mathcwson, Ford & Buchhorn 1992). In this reconstruction, backside infall can clearly be seen for both 
Centaurus/GA (centred at SGX ~ -3500km s"\ SGY ~ 1500km s _1 in figure 4c) and P-P (centred at 
SGX ~ 5500km s _1 , SGY ~ —1500km s _1 in figure 4c). A cone including the Centaurus/Great Attractor 
region (305° < I < 325°, 10° < b < 25°) is shown in §6 as figure 9c. This shows velocity turnaround at 
r ~ 3500 km s _1 , with line-of-sight infall velocity reaching a maximum of <~ 300 km s _1 at r ~ 4500 km s _1 . 

5. Effects of Parameters on Reconstruction 

In order to assess the effects of the parameters discussed in §3, the reconstruction was also carried out 
with a variety of different parameter values. When expanding the initial, redshift-space density field in terms 
of harmonics, the expansion at some radial and angular mode. This determines the resolution of the initial 
expansion, and hence the final reconstruction, as the whole procedure is carried out in terms of these same 
harmonics. If the expansion is terminated at a very low order, structures on a smaller scale than wavelength 
of the highest harmonic will be lost. On the other hand, if the expansion is continued to a very high order, 
we risk introducing artefacts as the density field ripples on scales smaller than the real structures. Appendix 
B of FLHLZ contains a detailed discussion of these problems. 

Figure 6 shows the density field on the SGZ = km s _1 plane, demonstrating how additional structure 
is resolved as angular modes are added. In figure 6a, the angular modes are terminated at Z max = 4, and 
the circular structures associated with low-frequency angular modes are very clear. Overdensities occur 
on distinct shells, and follow circular paths about the origin; furthermore, resolution falls with radius, as 
the wavelength of the highest order mode increases. In figure 6b, Z max = 10, and the structures are much 
more detailed, with less obvious artefacts. The arc across the top right still follows a circular path, and do 
the two clusters at the bottom left. Nonetheless, structure is resolved quite well, with both the P-P and 
Centaurus/GA regions clearly visible. Comparing these to the reconstructions earlier in this paper, where 
Imax — 15, it is interesting to note that the purely angular arc across the top left of the SGZ = km s _1 
slice is still present. 

Beyond those parameters inherent in the expansion, a choice of cosmological parameters is necessary for 
the actual reconstruction. The shape of the prior power spectrum, characterised by T and as, seems not to 
affect the reconstruction strongly within r < 10,000 km s . We have also performed the reconstruction of 
the IRAS velocities with a standard CDM (r = 0.5) prior. Since the standard CDM model has less large 
scale power than the T = 0.2 model, the WF smoothes more on large scales and therefore the reconstructed 
velocities tend to be smaller; the overall difference is however small with (DeltaRMs) ~ 50 km s _1 . Figure 
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7 shows density and velocity fields on the SGZ = km s -1 plane for two different prior power spectra; 
(r = 0.2, cr 8 = 1.0) and (T = 0.5, erg = 0.7). Although the choice of (3 plays a complex role in the mechanics 
of the reconstruction, changes within the range 0.2 < (3 < 1.0 have little effect other than a linear scaling of 
peculiar velocities. Figure 8 shows the density fields on the SGZ = km s _1 plane for (3 = 0.5 and (3 = 1.0, 
and the residual differences in the density and velocity fields. The residual velocity field is calculated by 
correcting for linear velocity scaling such that |v| oc (3. Hence: 

vi v 2 

Vresid = -7. ~T ■ {< ) 

Pi P2 

As can be seen from the plots, there is very little discrepancy from this linear scaling, indicating that within 
this volume, the reconstruction is robust to changes in (3. 



6. Comparison with Mark III 

The Mark III catalogue (Willick et al. 1995) contains about 3000 galaxy peculiar velocities from Tully- 
Fisher distances. In figure 9, we show comparisons between Mark III velocities and those reconstructed 
from the IRAS 1.2Jy survey. These are plotted for cones selected by Galactic latitude and longitude, 
following Faber & Burstcin (1988). Throughout these plots, it is important to note that the Mark III data 
is unsmoothed and hence displays considerably more scatter than that from WF reconstruction. As such, 
it is not possible to make anything more than a qualitative comparison. Davis, Nusser & Willick (1996) 
recently made a mode-by-mode comparison of IRAS and Mark III using spherical harmonics, finding a 
reduced \ 2 ~ 2. A summary of other comparisons is given in Dekel (1994). 

Figure 9a shows the central region of the Virgo cluster (265° < I < 315°, 67° < b < 80°). There is good 
agreement on the velocity gradient at around 1500 km s , where data are rich in both catalogues. However, 
for 2000 km s _1 < r < 3000 km s _1 , the WF velocities are near zero, while Mark III shows a continuation 
of the earlier velocity gradient, giving strong backside infall (v ra diai — —1200 km s _1 at r ~ 2500 km s _1 ). 
A similar effect is seen in figure 9b, which shows the Leo Cloud (200° < I < 260°, 50° < b < 70°). Velocity 
gradients in the two catalogues are very similar out to about 1500 km s _1 , after which the WF levels out 
while the Mark III continues to show back-side infall out to r ~ 3000 km s _1 ; beyond this, Mark III data is 
dominated by scatter. The Centaurus/Great Attractor region (305° < I < 325°, 10° < b < 25°) is shown in 
figure 9c. In both surveys clear back-side infall can be seen, with turn-around at r ~ 3500 km s _1 . Apart 
from the expected scatter in the Mark III, there is very good agreement. Figure 9d shows the Fornax- 
Eridanus region (193° < I < 245°, —66° < b < —46°). Again, the Mark III data shows considerably higher 
velocity contrast than the WF reconstruction, with a <~ 1000 km s _1 velocity difference across the nearby 
cluster, compared to <~ 300 km s _1 for WF. 

Many of these differences might well be inherent to the WF procedure, which, by construction, goes to 
Vpcc — at large scales. In regions of sparse IRAS data, this could account for the reduced velocity contrast 
in the WF reconstruction. In principle, a comparison between reconstructed velocities and those from direct 
measurement could be used to constrain power from larger scales. 



7. The Acceleration of the Local Group 



The dipole temperature anisotropy of the microwave background (CMB) has been measured by the COBE 
satellite to extraordinary accuracy, D = 3.343±0.016 mK in the direction (I = 264.4°±0.3°, b = 48.4°±0.5°) 
(Smoot et al. 1991, 1992; Kogut et al. 1993; Fixscn et al. 1994). Although various alternative theories have 
been proposed (e.g., Gunn 1988; Paczyhski & Piran 1990; Langlois & Piran 1996), this anisotropy is usually 
interpreted as due to the motion of the Earth with respect to a rest frame defined by the CMB. After 
correcting the Earth's motion relative to the Local Group (LG) barycenter, one infers (Smoot et al. 1991, 
Kogut et al. 1993) that the LG is moving at a velocity of 627 ± 22 km s" 1 in a direction (/ = 276° ±3°, b = 
30° ±3°). 
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In the gravitational instability scenario, the peculiar velocity of the LG is generated by surrounding 
fluctuations in the mass density. In linear theory, the relation between the gravitational acceleration induced 
by the mass inhomogeneities and the LG dipole is particularly simple, 



In the context of the linear biasing model adopted in our reconstruction procedure, the the integral over 
the mass fluctuations can be replaced by an integral over the galaxy fluctuations provided we allow for a 
relative bias, i.e., /(O) — > (3 = f(fl)/b. Note the dipole is independent of the extragalactic distance scale if 
we measure distances in km s _1 . 

The integral in equation 8 can be evaluated using a galaxy redshift survey. By comparing the result 
with the known LG velocity, one infers an estimate of (3. This technique (and variants using a flux-weighted 
dipole) has been applied by many authors to both optical and infrared galaxy samples as well as clusters 
catalogues (see the review by Strauss & Willick and references therein). The main uncertainty in the value 
of (3 determined in analyses of the dipole stem from statistical noise due to sparse sampling, the effects of 
redshift distortion, and the unknown contribution to the dipole from scales larger than the sample size or 
structures hidden behind the Zone of Avoidance. 

The dipole predicted from equation 8 is (with our prescription of linear bias) linear in the galaxy density 
field. Consequently, the dipole computed from our WF algorithm will be an optimal (in the sense of minimum 
variance) estimate of the dipole due to matter within the reconstruction volume to the extent that non-linear 
effects can be neglected. Figure 10 shows the predicted dipole generated by matter within radius R as a 
function of R as computed from our WF reconstruction (see Appendix A) for our canonical prior (r = 0.2, 
(Ts = 0.7, and (3 = 0.7). The amplitude of reconstructed dipole for the model is slightly larger than the 
dipole observed by COBE. The amplitude reconstructed dipole is roughly proportional to (3 (there is some 
non-linear dependence which enters in the correction for the redshift distortion, see equations 16 and 18 
of FLHLZ) and therefore a reconstruction with value of (3 slightly smaller than our canonical model would 
provide a better fit. The direction of the LG dipole is shown in figure 11. The misalignment angle between 
the IRAS and CMB Local Group dipoles is only 13° out to R = 50 /i _1 Mpc, partially due to the 'tug-of-war' 
between the Centaurus/GA and P-P superclusters. However, at R = 200 /i _1 Mpc, the dipole of our canonical 
reconstruction points in the direction (I = 247°, b = 37°) which is 25° away from the COBE direction. As 
our procedure suppresses shot-noise this misalignment may indicate non-linear effects on these very large 
scales. However, assessing the significance of this misalignment is difficult, as the dipole observed from any 
volume-limited sample will inherently show scatter about the 'true' CMB value (Vittorio & Juszkiewicz 1987, 
Juszkiewicz, Vittorio & Wyse 1990, Lahav, Kaiser & Hoffman 1989). For example, Lahav et al. (1989) show 
that for a biased CDM universe, 85% of observers will see a misalignment angle of 9 < 20° in a sample out 
to r = 4000 km s" 1 . 

In figure 11, wc show the scatter in the reconstructed dipole due statistical noise and from contributions 
from fluctuations outside our reconstruction volume. The statistical scatter about the dipole can be computed 
within the framework of the WF algorithm and the result for the scatter in the dipole is given in Appendix B. 
This scatter is shown in figure 1 1 as the set dotted curves; this scatter grows with distance as the smoothing 
by the WF increases to mitigate the effects of shot noise. The uncertainty in the dipole due to matter outside 
the reconstruction volume can be calculated given the power spectrum assumed in the prior; the result is 
given in Appendix B and shown (for our canonical prior) as the dashed lines in figure 11. This scatter 
decreases steadily to zero as R increases. We should point out that the dipole generated within r < R is a 
well defined quantity. 



8. Bulk Velocities 




(8) 
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A simple statistic of the peculiar velocity field is the average or bulk velocity within a window, W(r), 



A commonly used window is spherical top-hat, W(r) — (47r/3i? 3 ) _1 9(r — R), where 0(x) is the usual step 
function. The bulk velocity is sensitive to the fluctuations on scales > R and thus probes the linear regime 
even when non-linear effects in the density field on scales < R are important. 

It is straightforward to relate the bulk velocity of spherical window centred on the LG to the expansion 
of the density field in spherical harmonics (see Appendix C). Once again, the WF reconstructed velocity 
field can be used to compute an optimal estimate of the bulk velocity in linear theory. As in the dipole 
calculation there will be scatter introduced by dilute sampling and fluctuations outside the reconstruction 
which can be computed for a given prior model of power spectrum (see Appendix C). 

In figure 12, we show the Cartesian components of the bulk flow computed within spherical windows of 
radius, R. The points with error bars show the WF reconstruction. The curves show the values from the 
Potent analysis (Dekel 1994 and private communication) of the Mark III compilation of peculiar velocities 
(Willick et al. 1995). The reconstructed IRAS bulk flow out to 5000 km s _1 is ~ 300 km s _1 , which agrees 
in amplitude with that derived from the Mark III peculiar velocities (~ 370 km s _1 ). However, the two 
vectors deviate by ~ 70° in direction; there is a 33 percent probability of two random vectors aligning within 
70° . This discrepancy could be due to differences in the width of the zone of avoidance and other systematic 
effects. 

9. The Extent of the Supergalactic Plane 

The so-called Supergalactic Plane (SGP) was recognised by de Vaucouleurs (1956) using the Shapley- 
Ames catalogue. This followed earlier work by Vera Rubin, whose analysis of radial velocities of nearby 
galaxies suggested a differential rotation of the 'metagalaxy'. Indeed, the SGP had in fact already been 
noted by William Herschel more than 200 years earlier. Traditionally the Virgo cluster was regarded as the 
centre of the Supergalaxy, and this was termed the 'Local Supercluster'. However, recent maps (e.g. figures 
2-5) of the local universe indicate that much larger clusters (e.g. Great Attractor, Perseus-Pisces) are major 
components of this 'plane'. The north pole of the standard SGP (de Vaucouleurs 1976) lies in the direction 
of Galactic coordinates (I = 47.37°, b = 6.32°). The origin of SGL is at (I = 137.37°, b = 0°), which is one 
of the two regions where the SGP is crossed by the Galctic Plane. The Virgo cluster is at SGP coordinates 



Although the SGP is clearly visible in whole-sky galaxy catalogues, it has only been re-examined quan- 
titatively in recent years. Tully (1986) claimed that the flattened distribution of clusters extended across a 
diameter of <~ 0.1c with axial ratios of 4:2:1. Shaver & Pierre (1989) found that radio galaxies were more 
strongly concentrated towards the SGP than were optical galaxies, and that the SGP (as represented by 
radio galaxies) extended out to redshift z <~ 0.02. Di Nella & Paturel (1995) examined the SGP using a 
compilation of nearly 5700 galaxies larger than 1.6 arcmin, and found qualitatively agreement with the stan- 
dard SGP. Lahav et al. (1996, in preparation; also Lahav 1996) revisited the SGP using the Optical Redshift 
Survey (Santiago et al. 1995) and the IRAS 1.2 Jy survey (Fisher et al. 1995) To objectively identify a 'plane' 
they calculated the moment of inertia (with the observer located at the centre) by direct summation over 
the galaxies (taking into account the selection function) and subtracting the mean background density nb g 
in the absence of the 'plane'. By finding the eigenvalues and vectors of the inertia tensor they deduced that 
the derived 'plane' is aligned to within 30° of the standard SGP. The diameter of the SGP in ORS and IRAS 
was estimated to be at least 12,000 km s^ 1 . 

Here we take a different approach as the reconstructed field is continuous. Again, we can use our formalism 
to obtain the optimal reconstruction (in the minimum variance sense) of the moment of inertia. As explained 




(9) 



r<R, 



(SGL = 104°, SGB = -2°). 
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earlier, a convenient property of this Wiener approach applied to the density fluctuation field, S, is that it 
will also give the optimal reconstruction for any property which is linear in S. In particular, if we seek the 
optimal reconstruction of the moment of inertia 



Cij — Cij 



Ci 



1 

N 



[p (x) - n b g] {xi - Xi) (xj - Xj) dV 



This can be re-written as 



Cij — 



4ttR 3 



hi + 



where S k is the Kroneker delta and 



lij 



l 5(r)x l Xj 

JR 



P 
dV 



R 2 . 



(10) 



(11) 



(12) 



This can be expressed analytically in terms of the reconstructed coefficients 5f^ nn . The full mathematical 
details are given in Appendix D. We emphasise again that in the Wiener approach the density field goes to 
the mean density at large distances. This does not necessarily mean that the SGP itself disappears at large 
distances, it only reflects our ignorance of what exists out there, where only very poor data are available. 

To find the alignment and extent of the 'plane', we diagonalise the covariance matrix and find the eigen- 
values and eigen-vectors : 

Cu a = \ a u a , (13) 

where the A a 's and u a 's are the eigen-values and eigen-vectors respectively (a = 1,2,3). The 'half-width' 
(1-sigma) along each of the 3 axes is given by \Aq- Note that since the background contribution (the last 
term in equation 11) is isotropic, it only affects the eigen-values, but not the directions of the eigen-vectors. 

We have applied this technique to the reconstructed, real-space density field, assuming as priors = 0.7, 
a 8 = 0.7 and a CDM power spectrum with T = 0.2. Out to i? max = 4000,6000,8000 km s" 1 the derived 
plane is aligned with de Vaucouleurs' SGP to within 15°, 35° and 31° respectively. Note that the probability 
of 2 random vectors being aligned within an angle 8 is 



P(<6) 



1 - cos(0) 



(14) 



Hence, for 6 = 30° only P ~ 7%. 

Out to i? max = 4000 km s _1 the axial ratio is roughly 7:5:1, assuming a background ratio (ribg/p) = 
0.79, derived from the IRAS data. The results are only slightly different for raw harmonics, uncorrected for 
rcdshift distortion and noise, suggesting that redshift distortion is negligible on these large scales. The results 
also agree well with those direct summation of the moment of inertia (Lahav et al. 1996, in preparation). 



10. Conclusions & Discussion 



In this paper, we have applied Wiener reconstruction with spherical harmonics and Bessel functions to 
the IRAS 1.2Jy rcdshift survey. Using a prior based on parameters — 0.7, a 8 — 0.7 and a CDM power 
spectrum with T = 0.2, we find that: 

1. The reconstructed density field clearly shows many known structures, including N1600 (Santiago et 
al. 1995) and A3627 (Kraan-Korteweg et al. 1996). Notably, Perseus-Pisces is seen to extend out to 
R <~ 9000 km s _1 , while Virgo, Centaurus and Telescopium-Indus-Pavonis join to form a large structure 
extending over 1500 km s _1 < R < 7000 km s _1 . A number of new clusters are also observed. Detailed 
maps are shown in §4. 

2. The reconstructed velocity field shows backside infall for both Perseus-Pisces and the the Centau- 
rus/Great Attractor region. The Centaurus/Great Attractor supercluster shows velocity turnaround 
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at r ~ 3500 km s -1 , with line-of-sight infall velocity reaching a maximum of ~ 300 km s _1 at r ~ 
4500 km s _1 . Further, there is reasonable qualitative agreement between the reconstructed velocity field 
and that derived from Tully-Fisher measurements (Mark III). 

3. The misalignment angle between the CMB and reconstructed IRAS LG dipoles falls to a minimum of 
13°, calculated for R < 50 h^Mpc, but increases to 25° for R < 200 h^Mpc. 

4. The reconstructed IRAS bulk flow out to 5000 km s _1 is <~ 300 km s _1 , which agrees in amplitude with 
that derived from the Mark III peculiar velocities (~ 370 km s _1 ). However, the two bulk flow vectors 
deviate by some ~ 70°. 

5. The alignment and extent of the reconstructed Supergalactic Plane can be determined analytically from 
the harmonic coefficients. Out to R max = 4000, 6000, 8000 km s _1 the derived plane is aligned with de 
Vaucouleurs' SGP to within 15°, 35° and 31° respectively. Out to i? max = 4000 km s _1 the axial ratio 
is roughly 7:5:1. 

We also confirm that the reconstruction is robust to changes in prior parameters. 

Wiener reconstruction is particularly well suited to recovering real-space density and velocity fields from 
near whole-sky surveys. The WF formalism provides a rigorous methodology for variable smoothing, de- 
termined by the sparseness of data relative to the expected signal. As such, natural continuations of this 
work would be application to newer surveys, such as the PSCZ (Saunders et al. , in preparation), which is 
complete down to 0.6Jy, or ORS (Santiago et al. 1995). 
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Appendices 



Appendix A: Wiener Filtering 



Here, we give a brief review of the Wiener Filter technique; the reader is referred to LFHSZ, ZHFL, and 
Rybicki & Press (1992) for further details. Let us assume that we have a set of measurements, {d a } (a = 
1, 2, ... TV) which are a linear convolution of the true underlying signal, s a , plus a contribution from statistical 
noise, ep, such that 

d a = TZ a p [sp + ep] , (Al) 

where lZ a p is the response or "point spread" function (summation convention assumed). In the present 
context it would be the radial coupling matrix discussed in the previous section. Notice that we have assumed 
that the statistical noise is present in the underlying field and therefore is convolved by the response function. 

The WF is the linear combination of the observed data which is closest to the true signal in a minimum 
variance sense. More explicitly, the WF estimate is given by s a (WF) — F a p dp where the filter is chosen to 
minimise (\s a (WF) — s a \ 2 ). It is straightforward to show (see ZHFL for details) that the WF is given by 

F a p = (s a d 1 )(d 1 dg}- 1 , (A2) 



where 



(s a dl) = Tlp 1 (s a s\) 
(d a dl) = TZ ai Tips (s 7 4) + ( e 7 e I) 



(A3) 
(A4) 



In the above equations, we have assumed that the signal and noise are uncorrelated. From equation A4, it 
is clear that in order to implement the WF one must construct a prior which depends on the variance of the 
signal and noise. 

The dependence of the WF on the prior can be made clear by defining signal and noise matrices given 
by S a p — (s Q sJj) and N a p = (e a ejj). With this notation, we can rewrite equation A4 as 



s(WF) = S[S + N] _1 7e _1 d 



(A5) 



Formulated in this way, we see that the purpose of the WF is to attenuate the contribution of low signal 
to noise ratio data and therefore regularize the inversion of the response function. The derivation of the 
WF given above follows from the sole requirement of minimum variance and requires only a model for the 
variance of the signal and noise. The WF can also be derived using the laws of conditional probability if the 
underlying distribution functions for the signal and noise are assumed to be Gaussian; in this more restrictive 
case, the WF estimate is, in addition to being the minimum variance estimate, also both the maximum a 
posterior estimate and the mean field (see LFHSZ, ZHFL). For Gaussian fields, the mean WF field can be 
supplemented with a realisation of the expected scatter about the mean field to create a realisation of the 
field; this is the heart of the "constrained realisation" approach described in Hoffman & Ribak (1991; see 
also ZHFL). 

As Rybicki & Press (1992) point out, the WF is in general a biased estimator of the mean field unless the 
field has zero mean; this is not a problem here since we will perform the filtering on the density fluctuation 
field which has, by construction, zero mean. 



Appendix B: Dipole Velocity in Spherical Harmonics 
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The velocity of the central observer in the CMB frame, commonly referred to a "dipolc" , due to the 
fluctuations with r < R can be related to the spherical harmonic expansion of the density field (see FLHLZ 
equation C9), 

O - 6 r r' 

v(0) r< * = — Jd^'Siv')- (Bl) 

Vr 

R 

= ^ E^™^/ dr 'Mk n r') J dat'Y lm {v') 

Iran n Q ^ 

1 /2 

= £ (f ) E IT {-^ Re ^} * + V^/m[Jf ln ] y + J?e[<5f 0n ] z) (1 - j (k n R)) . 

^ ' n n 

Here i?e[a] and Im[a] refer to the real and imaginary parts of a complex number, a. The last line provides 
a convenient expression for the Cartesian components of the dipolc. In the limit that R — > oo, we recover 
the true dipolc, apart from statistical noise and systematic errors introduced by non-linear effects. ^From 
the last line equation (Bl), we see that the dipole within r < R can thus be regarded as the sum of two 
terms: the true CMB dipole and a correction for finite sample size (the jo(k n R) term). In fact the latter 
term involving the jo(k n R) can be shown to be the velocity of a shell at a distance R in the CMB frame. 

There are two main sources of error in the dipole derived via equation (Bl). First, there is a statistical 
noise introduced by the finite sampling of the density field by the IRAS catalogue. The scatter in the 
reconstructed dipole can be easily calculated in the framework of the Wiener filter (see FLHLZ Appendix F 
for the scatter in the reconstructed density and velocity field). The scatter in each component of the dipolc 
is given by 

(A* = ^ E *n C W [(I _ Fi) {m) 

nn' 

where F; = S/(S; + N;) _1 is defined in terms of the signal, S/, and noise, Nj matrices. 

The second source of error is from fluctuations outside the volume used in the Wiener reconstruction. 
The contribution to the dipole from r > R, 

vo«*(0) = £ / d?v'5(v')^ , (S3) 

r<R. 

introduces a 1-D rms scatter in the dipolc 

(Av(0))L = \^2 J dkP(k) [ ]o (kR)f . (54) 

In the above equation, P(k) is the linear power spectrum defined by the prior. The total error in the 
reconstructed dipole is taken to be the quadrature sum of equations B3 and B4. 



Appendix C: Bulk Velocity in Spherical Harmonics 



The average or "bulk" velocity measured within a spherically symmetric region defined by window func- 
tion with characteristic scale r = R s is defined by 

(v) fl . = J d 3 rW(r)v(r) , (CI) 
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where the volume integral of the window is taken to be unity. If we rewrite the above expression in terms of 
the Fourier components of v, 

W) Rs = j^ J d 3 kv k W(kR s ) , 

where W(kR s ) is the Fourier transform of the window function (W(kR s ) = 3ji(kR s )/kR s for spherical 
top-hat). From the above equation, we see that the bulk velocity for region centred on the LG is formally 
equivalent to the dipole velocity of LG induced by the density field obtained from smoothing with the 
window. Consequently, the analogous formulae for the bulk velocity, uncertainty due to sampling, and the 
contributions from fluctuations outside the reconstruction volume in terms of the spherical harmonics are 
readily derived by substituting 

Simn^Si mn W(k n R s ) and P(k) -» P(k)W(k n R s ) 

in equations Bl, B2, and BA 



Appendix D: Moment of Inertia in Spherical Harmonics 



Consider the moment of inertia tensor within a sphere of radius R around the origin: 

la = I S(r)xiXjdV . (1) 



hj — I <■' \ * /•'/■<; 1 
JR 



We now express the density fluctuation in terms of spherical harmonics and Bessel functions: 

Z max m=+l n„,x(I) 



I m= — l n=l 



By substituting eq. (2) into eq. (1) and utilising the properties of the spherical harmonics we find after some 
algebra the elements of the symmetric matrix 7jj : 



"max (2) 



n=l 



(#2,-2,71 + °~2 : 2,n) — \ ~7= S 2 fi.n 



n»ax(0) 



^ Co,ri-Bo,n<5o,0,r 



n=l 



"max (2) 

lyy = £2,71-82,71 



n=l 



/2tt , s Air 

(02, -2, n + 02,2,n / A/ 45 2 '°' r ' 



, ^ n ma *(0) 

+ —Z— ^ Co,rt-So,n^O,0,n 
ra=l 



<(2) 



'nmax(O) 



Co,n-Bo,n^O,0,-, 



n=l 



-«max(2) 



— IA 



^ C2, n B2,n{^2,-2,n — $2,2,n) 



71=1 



k(2) 



C2,nB2,n{^2,-l,n ~ 0~2,l,n) 



n=l 



lyz — C2,n^2,n(^2,-l,n + °~2,l,ra) 



c(2) 



71=1 
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where n max (/) is the n max corresponding to a particular value of I, and 

f R 

Bin = / ji(k n r)r A dr 
Jo 



Note that as the Uj is quadratic in the coordinates, the only harmonics to appear are of I = and 1 = 2. 
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Figure Captions 



Fig. 1 — Scatter between velocities generated in Local Group and CMB frames, plotted as a function 
of used in the reconstruction. The amplitude of the scatter is divided by so as to remove the linear 
scaling of velocities with inherent in reconstructions based on linear theory. A clear minimum can be seen 
at ~ 0.7. 

Fig. 2 — The reconstructed density field, evaluated on thin shells at various real space distances, shown 
in Galactic Aitoff projection. Dashed lines show 5 < 0, and solid lines show 5 > 0, with contour spacing 
AS = 0.1. (a) 2000 km s -1 : Clusters marked are N5864, Virgo (Vir), Centaurus (Ccn), Hydra (Hyd), 
Fornax- Doradus-Eridanus (FDE), Ursa Major (Urs), and Ca. Voids marked are Local Void (LV), Va 
and Y0. (b) 4000 km s _1 : Clusters marked are Centaurus (Cen), Telescopium-Indus-Pavonis (TIP), Hydra 
(Hyd), Leo, Cancer (Can), N1600, Camelopardalis (Cam), Perseus-Pisces (P-P), Cetus (Cet), C0, C7 and 
CS. Voids shown are Local Void (LV) and Y0. (c) 6000 km s" 1 : A3627, Hydra (Hyd), Leo, Cancer (Can), 
Orion (Ori), A569, Perseus-Pisces (P-P), Pegasus (Peg), Telescopium-Indus-Pavonis (TIP), C7, CS, Ce and 
CC One void is marked V7. (d) 8000 km s" 1 : Clusters marked are Great Wall (GW), A3627, Hydra (Hyd), 
Leo, A539, A779, A400, Coma (Com), Perseus-Pisces (P-P), Pegasus (Peg), Cygnus (Cyg), C7, CC and Cr). 
The continuation of V7 is shown. 

Fig. 3 — The reconstructed velocity field, evaluated on thin shells at various real space distances, shown in 
Galactic Aitoff projection. These velocity fields correspond with the density fields shown in figure 2. Dashed 
lines show infall, and solid lines show outflow. The first solid line is for u ra diai = km s _1 , and contour 
spacing is Av ra di a i = 50 km s _1 . (a) 2000 km s _1 : Note the velocity dipole due to the Centaurus/GA region 
(towards I ~ 300°, b ~ 15°). (b) 4000 km s" 1 : Note outflow towards P-P (I ~ 145°, b ~ -20°), N1600 
(/ ~ 195°, b ~ -20°) and Cancer (I ~ 190°, b ~ 20°). (c) 6000 km s _1 : Dominated by motion from the 
south Galactic pole towards the North Galactic pole, caused by outflow towards the Great Wall in the 
northern Galactic hemisphere, (d) 8000 km s _1 : Outflow caused by large superclusters at I ~ 60°, b ~ 60° 
(forming behind the Great Wall), and I ~ 220°, b ~ -30° (forming behind Cr) and A400). Also, note strong 
infall can be seen behind Pegasus (I ~ 90°, b ~ -25°) and Hydra (/ ~ 285°, b ~ 5°). 

Fig. 4 — The reconstructed density field on thin slices at Supergalactic x,y,z — 0, ±4000 km s _1 . 
Contours start at S = with spacing Ad — 0.25. Note the clear presence of the Supergalactic plane in the 
SGZ = km s" 1 plot (e). 

Fig. 5 — The reconstructed velocity field on thin slices at Supergalactic x,y,z = 0, ±4000 km s _1 . 
Arrows show the projection of the local peculiar velocity onto the plane, with length giving the amplitude in 
accordance with the axis scale. These velocity fields correspond to the density fields shown in figure 4. Note 
clear backside infall towards both Centaurus/GA (centred at SGX ~ —3500 km s _1 , SGY ~ 1500 km s _1 
in (e)) and P-P (centred at SGX ~ 5500 km s"\ SGY ~ -1500 km s" 1 in (e)). 

Fig. 6 — Cone diagrams comparing the IRAS reconstructed velocity field with that derived from Tully- 
Fisher measurements (Mark III). Hollow triangles show the IRAS reconstruction, while filled squares show 
Mark III. (a) Central region of the Virgo cluster (265° < I < 315°, 67° < b < 80°). (b) Leo Cloud 
(200° < I < 260°, 50° < b < 70°). (c) Centaurus/Great Attractor region (305° < I < 325°, 10° < b < 25°). 
(d) Fornax-Eridanus region (193° < I < 245°, -66° < b < -46°). Note that the Mark III data is unsmoothed 
and hence displays considerably more scatter than that from WF reconstruction. 

Fig. 7 — Demonstration of different reconstruction resolutions. The reconstructed density field, evaluated 
on a thin slice at SGZ = km s _1 . Contours start at 6 = with spacing AS = 0.1. Figure 5c shows the 
corresponding plot for l max = 15; the default used throughout this paper, (a) Z max = 4: Clear circular 
patterns can be seen in the density field, (b) Z max = 10: Only a few circular artefacts remain, given the 
addition of higher-order harmonics. 

Fig. 8 — Effects of varying cosmological parameters T and as- (a) and (c) show the reconstructed density 
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field, evaluated for a thin slice at SGZ = km s _1 . Contours start at S = with spacing AS = 0.1. (b) and 
(d) show the reconstructed velocity field, evaluated for the same thin slice. Arrows show the projection of 
the local peculiar velocity onto the plane, with length giving the amplitude in accordance with the axis scale. 
Figures 5e and 6e show the corresponding plots for = 0.7, T = 0.2, erg = 0.7; the canonical parameters 
used throughout this paper, (a) and (c) = 0.7, T = 0.5, cr 8 = 0.7. (b) and (d) = 0.7, T = 0.2, cr 8 = 1.0. 

Fig. 9 — Effects of varying the cosmological parameter 0. (a), (b) and (c) are density fields, evaluated 
for a thin slice at SGZ = km s _1 . Contours start at 5 — with spacing AS = 0.1. (a) and (b) show 
reconstructed density fields for of 0.5 and 1.0 respectively, with T = 0.2 and erg = 0.7. Figures 5e shows 
the corresponding plot for = 0.7, T = 0.2, erg = 0.7; the canonical parameters used throughout this paper, 
(c) shows the residual density field after subtraction of (a) from (b), demonstrating the extremely small 
differences caused by change of 0. Finally, (d) shows the residual velocity field, evaluated for the same 
thin slice. Arrows show the projection of the local peculiar velocity onto the plane, with length giving the 
amplitude in accordance with the axis scale. 

Fig. 10 — The amplitude of the acceleration, or dipole motion, of the LG caused by the fluctuations within 
a sphere of radius R versus R. The heavy solid curve shows the dipole derived from the Wiener reconstruction 
with our canonical prior (r = 0.2, erg = 0.7, and = 0.7). The dotted curves show the expected scatter from 
finite sampling (see Appendix B) while the dashed curves show the scatter due to fluctuations outside the 
reconstruction volume (i.e., r < 200 /i _1 Mpc); the dot-dashed curve is the quadrature sum of both terms. 
The dotted horizontal line is the value of the LG dipole inferred by COBE, 627 km s _1 . 

Fig. 11 — The direction of the LG dipole. The crosses show the convergence of the direction of the 
reconstructed dipole; starting at the top of the plot the crosses give the direction (in Galactic coordinates) of 
the dipole within R (see Fig. 1) as R is increased in 1 /i _1 Mpcintervals. The direction of the dipole inferred 
from the COBE measurements, (I = 276°, b = 30°; Smoot et al. 1991) is denoted by the star. The circular 
curves denote angular separations from the COBE result in 10° increments. 

Fig. 12 — The average or bulk velocity within a sphere of radius R centered on the Local Group a 
function of R. The four panels show the three Cartesian components and the scalar amplitude of the bulk 
flow. The points with error bars represent the bulk flows in the Wiener reconstruction with our canonical 
prior (r = 0.2, erg = 0.7, and = 0.7). The error bars are the scatter due to both the finite sampling 
and fluctuations outside the reconstruction volume (see Appendix C). The triangles connected by the curves 
represent the measurements from the POTENT reconstruction algorithm applied the Mark III peculiar 
velocity data (taken from Dekel 1994). The two curves represent two different weighting schemes and reflect 
the systematic uncertainty; the estimated random error is approximately 15% (Dekel 1994). 
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